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Abstract. In materials without spatial inversion symmetry the spin degeneracy of 
the conduction electrons can be lifted by an antisymmetric spin-orbit coupling. We 
discuss the influence of this spin-orbit coupling on the spin susceptibility of such 
superconductors, with a particular emphasis on the recently discovered heavy Fermion 
superconductor CePtsSi. We find that, for this compound (with tetragonal crystal 
symmetry) irrespective of the pairing symmetry, the stable superconducting phases 
would give a very weak change of the spin susceptibility for fields along the c-axis 
and an intermediate reduction for fields in the basal plane. We also comment on the 
consequences for the paramagnetic limiting in this material. 
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1. Introduction 

Our understanding of conventional superconductors is based on the BCS theory of 
Cooper pairing induced by electron-phonon interaction. The electrons pair up in the 
superconducting phase with the most symmetric pair wave function possible, i.e. in 
the s-wave orbital and spin-singlet channel. According to Anderson (1959) this type 
of pairing requires time reversal symmetry, because the paired electrons need to be in 
degenerate states of opposite momentum and opposite spin pQ. This theorem guarantees 
that the conventional superconducting phase remains stable even, if a sample is dirty as 
long as time reversal symmetry is present and allows us to construct degenerate electron 
pairs with vanishing total momentum. 

The situation becomes more complex for unconventional superconductivity with 
Cooper pairs of lower symmetry. Impurities of all kind are detrimental to such pairing 
due to the momentum averaging effect of random potential scattering. For clean samples 
new symmetry criteria have been formulated by Anderson (1984) |2j. While spin-singlet 
pairing only requires time reversal symmetry, spin-triplet pairing also requires inversion 
symmetry. These symmetries are present in most materials and no discussion on this 
issue is usually required. 

Motivated by the recent discovery of the heavy Fermion superconductor CePtsSi 
which has no inversion center, we will discuss the magnetic properties of the 
superconducting phase J3] While the absence of an inversion center does not imply 
immediately unusual behavior of the superconducting phase special interest in this 
case has arisen with the observation that the upper critical field dramatically exceeds 
the paramagnetic limit. For heavy Fermion superconductors the coherence length is 
generally rather short, roughly 100 A, giving rise to large orbital depairing fields. Thus, 
if there were a paramagnetic limiting effect, then it would be likely be observed; as for 
example in CeCoIns j3]. CePtsSi has a critical temperature of T c = 0.75K and a zero- 
temperature extrapolated upper critical field of 5 Tesla. The estimated paramagnetic 
limit is a factor 5 smaller than H c2 (H p ~ A/v^/Ub ~ 1 Tesla) [3]. Paramagnetic 
limiting is expected for spin-singlet pairing, as the Zeeman coupling to the spins would 
break up the Cooper pairs. Since Cooper pairs in spin-triplet configurations would 
not be destroyed by spin polarization, it might be concluded that spin-triplet pairing 
is realized in CePtaSi. However, the absence of an inversion center appears to be an 
obstacle for spin-triplet pairing. 

It has been shown that broken inversion symmetry is not indiscriminately 
destructive for the spin-triplet pairing, and, additionally, it softens the effect of the 
paramagnetic depairing for spin-singlet pairs p]. For many trpes of crystal lattices, 
the absence of an inversion center introduces an antisymmetric spin-orbit coupling 
(SOC), analogous to the well-known Rashba type of spin-orbit coupling. Naturally 
this influences the Cooper pairing through the modification of the band structure as 
shown by various groups EE Ej . 

We would like to discuss here the problem of the spin susceptibility of the 
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superconducting phase, since this is also directly connected with the issue of 
paramagnetic limiting. The spin susceptibility measured by means of the NMR-Knight 
shift in the superconducting phase is often used to distinguish between spin-singlet and 
spin-triplet pairing. We will show here that this kind of discrimination between the 
two types of states is not generally possible anymore. Rigorously speaking it is, of 
course, not possible to separate spin-singlet and triplet in the presence of SOC. We will 
consider here weak SOC which we can turn on adiabatically to follow the evolution of 
the originally well-defined singlet and triplet pairing states. Our basic result is that 
the spin susceptibility of the spin-singlet states gradually approaches the behavior of 
the spin-triplet state which survives the presence of the spin-orbit coupling. Moreover 
we can predict that the spin susceptibility for fields along the z-axis of the tetragonal 
crystal lattice of CePtsSi would be less suppressed than for fields along the basal plane. 



2. The basic Model and normal state Greens function 

For the following discussion we will use the model introduced by Frigeri et al. jHj for 
CePtaSi, which has the following single-particle Hamiltonian: 

n o=J2 t^ a ° + ag k ■ CT L' c k s c ks> (!) 
k,s,s' 

where ( c fcJ creates (annihilates) an electron with momentum k and spin s. The 
band energy £^ = — \i is measured relative to the chemical potential /i and ag k ■ er 
introduces the antisymmetric spin-orbit coupling with a as a coupling constant (we set 
(9%)k = 1 where () denotes the average over the Fermi surface). 

We give here a brief discussion on the origin and some basic properties of the 
antisymmetric SOC. In a crystal lattice the electrons move in a periodic potential U(r). 
In the absence of inversion symmetry there is no symmetry point in the unit cell relative 
to which U(r) = U(—r) is satisfied. This also implies that the Bloch function does not 
have the property that uy,{r) = u_ k (—r). Ignoring for the moment relativistic effects, 
the potential yields the following contribution to the single-particle Hamiltonian: 

i Ju.c. 

k,s 

= £/ <Pru_ k (r)U(r)u k (r)cl s c ks = Y, tj ^ k ) c i s c ks • ( 2 ) 

fcJ U.C. 7 

The integral runs over the unit cell of the lattice. The resulting potential U(k) is 
even in k, i.e. U(k) = U(—k). On this (non-relativistic) level the lack of inversion 
symmetry does not affect the band structure which remains symmetric under the 
operation k — ► —k due to time reversal symmetry, i.e. the fact that u,(r) = u_ k (r). 

Now we include the SOC. The symmetric SOC of each ion couples different atomic 
orbitals and requires a multi-orbital description. In this way the spins would be 
converted into pseudospins which can be handled formally in the same way as the original 
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spins. For our purpose the most important influence of this would be an anisotropic 
g-tensor which could be straightforwardly introduced into the following discussion and 
would affect the susceptibility of the normal and superconducting phase in the same way. 
For the sake of simplicity, however, we will keep here the g-tensor isotropic. Moreover 
we ignore the effect of spin-orbit coupling on the pairing interaction. 

For a lattice without inversion symmetry SOC already appears on the level of a 
single-band model yielding the second term in Eq.(JTJ [8]. The vector function ag k is 
derived from the relativistic correction ^3 \ v x ^rU(r)] ■ S, which yields 

a9k = ~2^J UC dh{J k^ x VrU(r)} (3) 

with 

Jfc(r) = \- U k( r ^ tk + V r)u k (r) + u k (r)(ik - V r )u%(r)] . (4) 

It is easy to verify that ag k = 0, if U(r) = U(—r) and u k (—r) = u_ k {r). In the 
absence of inversion symmetry, however, g k is finite and satisfies g k = —g_ k , since 
J-k( r ) = ~Jk( r ) ( n °te if g k = g_ k this implies time reversal symmetry is broken). 
With these properties it is now clear that the Hamiltonian Tio is invariant under time 
reversal T but not under inversion operation X, because 



l{ag k ■ <x}X 1 = —ag k ■ er and T{ag k ■ cr}T 1 = ag k ■ cr. 

The SOC yields a modified band structure. We parameterize the Green's function 



by 



G (k,iu n ) = G° + (k,iuj n )a + {g k ■ &)G Q _{k,iuj n ) 



where 



G±(k,iuj n ) 



± 



iuj. 



(5) 



(6) 



£ k -a\g k \ iu n -£ k + a\g k 

and g k = g k /\g k \ (\g\ = V» )■ The band splits into two spin dependent parts with 
energies E k± = £ k ± a\g k \. The spinor is twisted on the two bands in a way that is 
described by the antisymmetric part of the Green's function, (g k ■ &)G°_(k,iu n ). 



3. Superconducting phase 

Now we turn to the superconducting phase and introduce the general pairing interaction 

H pair = - £ V ^---^kM)cl J_ ks2 C_ kls3 C k ^ si (7) 
k , k , s 1 

The interaction satisfies the relations V^^fc, k') = — V^ a - tl> §(— k, k') = 
— Va,p;6,j{ki ~k') = Vs^-/3 >a (k', k) to assure that the Fermion sign and the time reversal 
symmetry (T) are preserved. 
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For the following calculations it is advantageous to discuss the superconducting 
phase by means of Green's functions 

GxJk, t) = 



~ Wc fc ^(0)}) 
( T Ac k ,x(T)c_ k ^0)}) 

where the operators A (r) and (r) are expressed in the Heisenberg representation. 
These Green's functions have to satisfy Gor'kov equations of the following form 

Go 1 (k,iu n ) + A(k)G^(-k,-iuj n )A^(k) G(k,iuj n ) = a . (8) 

F(k, iuj n ) = G (k, iu n )A{k)G T (-k, -iu n ). (9) 
F\k, iu n ) = G^(-k, -iuj n )A\k)G(k, iu n ) (10) 

where At (At) is the gap function defining the order-parameter of the superconducting 
state and is 2 x 2-matrix in spin space. 

In general the gap function has a singlet ip(k) and a triplet component d(k), i.e. 
At (As) = i{tp(k) + d(k) ■ (r}a y , where ip(k) = ip(—k) is even and d(k) = —d(—k) is odd. 
For a finite a, the gap equations for the spin-singlet and triplet channel are coupled 
[UJ Ej and thus the gap function is a mixture of both channels. However, this coupling 
is of the order a/ep which we take to be small (ep\ the Fermi energy, or analogue to 
the band width). Consequently, we ignore this coupling and consider the singlet and 
triplet channels separately. In particular, it manifests itself in two ways. It gives rise 
to different magnitude gaps on the two SOC spilt Fermi surface sheets. It should also 
be observable through a non-vanishing spin-current contribution in the excess current 
associated with Andreev scattering. 

It has been shown that once a » k^Tc for the spin-triplet channel, then the 
only spin-triplet state that is permitted satisfies d(At)||gr(At) 0. We assume that 
ksT c < a « €f and solve in sequence the equations (JB1 El HDj) . The solutions are 
formally the same for both singlet (A (At) = i{ip(k)}a y ) and triplet [A(At) = i{d(k)-&}a y 
with <Z(Ac)||gf(At)] pairing. We find 



G(k, iu n 
G±(iu n ) = 

where E± = £ ± a|<jr|, and 

F(k,iu n ) = 
P\k,iu n ) = 

F±(k,iu n ) = 



G + (k,u n )cr + (g k ■ &)G-(k,iu r , 



iuj 



± 



iu + E_ 



'u;2 + |A|2 + £2 



(11) 



[F + (k,iLU n )a + (9k ■ &)F-(k,iu n )]A(k), 
A\k)[F + (k,iu n )a + (g k ■ ar)F_{k,iu n )} 



± 



(ul + W + El) ^ M + IAI 



(12) 



Note that the anomalous Green's function has less symmetry than the gap function. In 
particular, A(At) is either symmetric (singlet) or antisymmetric (triplet) with respect to 
k — > —At, while the resulting anomalous Green's function has both a symmetric and an 
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antisymmetric component. This leads to the mixing of the spin-singlet and spin-triplet 
Cooper pairs and has been discussed by Gor'kov and Rashba in the case of 2D metals 
[Oj. However F_ is an odd function in £ and doesn't contribute to the gap equation 
in its weak coupling formulation. In this way the magnitude of the gap if) = ip(k) for 
the singlet s-wave order parameter can be approximated by the standard and universal 
BCS gap equation, 



ln ^ = " / dx rrr—2 1 r~r-^ (13) 

exp It k B T 

where T is expressed in units of T c , ip in units oii/j(T = 0), where tp(T = 0)/ ksT c = 
7r/7 and C = = 0.577 correspond to the Euler constant. The deviation will be 

of the order (a/e c ) 2 , with e c the cut-off energy of the attractive interaction, which 
is small in our case. For the protected spin-triplet state (d(k) = A g(k)) we have 
A (T = 0)/A;bT c ~ ip(T = 0)//cbT c , and to avoid numerical complication superfluous 
for the following discussion, we have set A (T/T C ) = ip(T/T c ). 



4. The Static Uniform Spin Susceptibility 



In materials with a spatial inversion center, the measurement of the Knight shift in 
the resonance frequency (5u) by nuclear magnetic resonance (NMR) is an important 
experimental tool in determining the nature of the superconducting state. In particular, 
it allows the determination of the spin structure of the Cooper pairs. The measurement 
of the temperature dependence of the ratio 5u s /5uj n = Xs/Xn is a direct measure of the 
behavior of the spin susceptibility of the superconducting state (xs) (Xn is the spin- 
susceptibility of the normal state). For spin-singlet superconductors it is known that 
the paramagnetic susceptibility is proportional to the density of normal electrons, which 
vanishes at zero temperature. In the spin-triplet case, the spin of the Cooper pairs can 
contribute to the susceptibility. In particular, if the external field (H) is parallel to 
the spin of the Cooper pair (H _L d), then the susceptibility coincides with that of the 
Fermi-liquid normal state. This property was used to experimentally confirm spin-triplet 
superconductivity, for example, in Sr 2 Ru0 4 [TU] . 

In principle, for materials with strong spin-orbit coupling, the total magnetic 
susceptibility cannot be split into separate orbital and spin parts. However, if a « ep 
it is possible to isolate the two components JT]. As shown in Ref. ^2], the spin 
susceptibility tensor xh i n the superconducting state can be expressed as 



Xtj = - /i|fci?T^^tr{6- i G(fc,^„)a i G(fc,u; ri ) - a, t F(k, u n )cr] P\k, u n )} 

k 

(14) 
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If we assume a spherical Fermi surface and a constant density of states N(£) close to 
the Fermi surface, then we can replace the sum over k by Ylk — > N(0) J ^ f d£. In the 
normal state the integral over d^ cannot be carried out before the sum over uj n because 
the regular Green's function will be formally divergent [TB] . Doing the sum first, we 
find that the normal state spin susceptibility corresponds with the Pauli susceptibility 
Xn — 2/x#iV(0). This would not be the case if electron- hole asymmetry is taken into 
account (i.e if dN(^)/d^=o ^ 0). For our purposes this is an unnecessary complication 
which we will neglect here. 

To avoid carrying out the summation over uo n in the superconducting state, we follow 
Abrikosov et al. and sum and subtract the expression corresponding to the normal state 
|13j . The integral of the difference between the integrands rapidly converge in this case. 
Consequently, the order of summation and integration can be interchanged. For a singlet 
gap function we find a generalization of the result obtained by Gork'ov and Rashba [E] 
and Bulaevskii et al [TT)] : 



X i 




I / i - g%,i \*P(k)\ 2 



+9k H" 2 n +\m\ 2 ) 3/2 /k 



(15) 



a 



For the triplet gap function (with d(k)\\gu) the susceptibility is independent of 

^-{'-^(^odS^U- (16) 

More precisely, the contribution due to a from the regular Green's function is cancelled 
out by the contribution of the anomalous Green's function. 

We now apply the results to the recently discovered heavy fermion superconductor 
CePtaSi 0. In this case, the generating point group symmetry is C& v for which the 
simplest form of is oc k x z = (—k y ,k x ,0) [3]. This has the same form as the 
well-known Rashba spin-orbit coupling [14] . The spin-susceptibility for the singlet s- 
wave gap function is shown in FigO The left plot shows the corresponding behavior 
of the susceptibility for the field along the c-axis (x\\ = Xc,c)- The right plot shows 
the spin susceptibility for the field in the ab-plane (x± = Xa,a = Xb,b) as function of 
the temperature for three different values of the spin-orbit coupling (a). In Fig. Owe 
show how the zero-temperature value of the susceptibility rises as a function of a. For 
a ^> k B T c the x\\ approaches Xn and x± — Xn/2. 

We remark that the susceptibility increases with the spin-orbit coupling strength. 
For a very large, the resulting susceptibility looks very similar to that obtained for the 
triplet p-wave gap function shown in Fig. El For the spin triplet phase we chose in 
accordance to Ref.|3] the pairing state d = xk y — yk x . The similar properties of the 
spin susceptibilities make it difficult to distinguish between a spin-triplet and spin-singlet 
order parameter through NMR measurements in the strong SOC limit. 

Due the complicated band structure of CePt 3 Si and the coexistence of 
superconductivity with antiferromagnetism [E], our theory does not provide a 
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Figure 1. The spin susceptibility in case of singlet s-wave gap function for 
9k 00 ("^ji 0) (CePtsSi). The spin susceptibility in the ab-plane x± an< i along 
the c-axis Xll as a function of T for three different values of the spin orbit coupling 
a. The susceptibility in the superconducting state (T/T c < 1) increases with the spin- 
orbit coupling strength. The susceptibility is more strongly suppressed in the ab-plane 
than along the c-axis. At T — we have x± — Xjj/2- 



x/x N 




Figure 2. The zero-temperature value of xjj/xn of the spin singlet state as a function 
of a. Note that xl = Xij/2. 



quantitative description for the spin susceptibility for CePtsSi. However, our approach 
illustrates the behavior expected at a qualitative level. More precisely, the susceptibility, 
independent of the kind of pairing and of the strength of the SOC, is more strongly 
suppressed in the ab-plane than along the c-axis. This angle dependence of the 
static uniform susceptibility should be confirmed by NMR-Knight-shift measurements. 
Moreover, our discussion supports the conclusion of Frigeri et al. that the spin "singlet" 
pairing state acquires a certain robustness against pair breaking due to spin polarization 
A rough estimate of the zero-temperature limiting field is obtained by comparison 
of superconducting condensation and spin polarization energy, leading to 

H p » kBTc . (17) 

/i BV /l-X s (T = 0)/ X n 

In principle this can become very big for fields along the c-axis and roughly 1- 2 Tesla 
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Figure 3. The spin susceptibility for a spin-triplet p-wave gap function d(k) \\ oc 
(— k y , k x , 0) (CePtsSi). The spin susceptibility in the ab-plane x± an d along the c-axis 
X|| as function of T. The susceptibility is in this case independent of the spin orbit 
coupling a. In the superconducting state the susceptibility in the ab-plane coincides 
with that of the normal state. 

for fields in the basal plane. 
5. Conclusions 

We have determined the spin susceptibility in superconductors without inversion 
symmetry. While the spin-triplet and spin-singlet order parameters are mixed in general, 
we can discuss predominantly spin-triplet or spin-singlet gaps when the spin orbit 
coupling is much smaller than the band width. We have found that for the surviving 
predominantly spin-triplet gap, the lack of inversion symmetry does not change the spin- 
susceptibility. For a predominantly spin-singlet gap, the lack of inversion symmetry leads 
to an increase in the spin-susceptibility. For large spin-orbit coupling (relative to T c ), 
the spin susceptibilities for both the spin-singlet and spin-triplet gaps become similar. 
For the heavy fermion superconductor CePt 3 Si, we have predicted that, independent 
of the pairing symmetry, the susceptibility is more strongly suppressed in the ab-plane 
than along the c-axis. More generally, we may state that the spin susceptibility in a 
superconductor without inversion center is approximately described by the behavior of 
the spin triplet superconductor with d(k) = A g^ and is obtained from Eq. (fTtjj) . 
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